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ABSTRACT 



Context. Recent X-ray observations have proved to be very effective in detecting previously unknown supernova remnant 
shells around pulsar wind nebulae (PWNe), and in these cases the characteristics of the shell provide further clues on 
the evolutionary stage of the embedded PWN. However, it is not clear why some PWNe are still "naked". 
Aims. We carried out an X-ray observational campaign targeted at the PWN G54.H-0.3, the "close cousin" of the Crab, 
with the aim to detect the associated SNR shell. 

Methods. We analyzed an XMM-Newton and Suzaku observations of G54. 1+0.3 and we model out the contribution of 
dust scattering halo. 

Results. We detected an intrinsic faint diffuse X-ray emission surrounding the PWN up to ~ 6' (~ 10 pc) from the 
pulsar, characterized by a hard spectrum, which can be modeled either with a power-law (7 = 2.9) or with a thermal 
plasma model {kT = 2.0 keV.) 

Conclusions. If the shell is thermal, we derive an explosion energy E = 0.5 — 1.6 x lO''^ erg, a pre-shock ISM density 
of 0.2 cm~^ and an age of ~ 2000 yr. Using these results in the MHD model of PWN-SNR evolution, we obtain an 
excellent agreement between the predicted and observed location of the shell and PWN shock. 

Key words. ISM: supernova remnants - (ISM:) dust, extinction. X-rays - ISM, X-rays: individuals: G54. 1-1-0.3 
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One of the most intriguing problems in the field of the 
Pulsar Wind Nebulae (PWNe) study is the lack of a 
shell around some of these objects. This is somehow dis- 
turbing for the consolidated picture of a remnant of a 
core-collapse supernova, which indicates that the PWN 
is expanding inside the host supernova remnant, giving 
rise to a variety of complex phenomena, like reverbera- 
tion, Ray leigh- Taylor instability at the interface between 
the PWN and ejecta, rejuvenating the shell, etc. (e.g. 
van der Swaluw et al.ll200lt iBlondin et al.ll200ll IChevalied 



2005t iGelfand et al.l l2009l and references therein). One of 



the reasons could be the lack of deep observations aimed at 
the PWN surroundings. Indeed, recently a shell-like com- 
pone nt has been observed in many objects, such as G21.5 - 



0.9 (iBandiera fc Bocchinol 120041: iBocchino et"all I2005D 



GO.9+0.1 (iPorauet et airi2003[) . 3C58 (|Bocchino et alJ 



I2OOII: iGotthelf et al.l 120071) . Therefore, X-ray observations 
are very effective for the discovery of associated shell com- 
ponents, even in the presence of high absorption column 
densities (G21. 5-0.9 has Nh ^ 2 x lO^^cm-^; GO.9+0.1 

-2\ 



even 



10^3 cm~^). 



The objects in which pulsar, plerion, and shell are all 
detected (collectively known as composite SNRs) are ex- 
tremely important to set the physical conditions for their 
modeling. The properties and the evolution of a PWN are 
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determined by the interaction of the pulsar wind with the 
ambient medium. The effectiveness by which the surround- 
ing matter confines the PWN is very important to deter- 
mine the the level of the synchrotron emission from the 
nebula. Therefore, measuring density and pressure in the 
shell component is needed for better constraining the mod- 
els of the PWN. Moreover, the pulsar, the PWN and the 
shell would allow us to estimate in independent ways some 
quantities, such as the actual age of the object, or the in- 
ternal pressure of the nebula. This redundancy would allow 
us also to verify our assumptions, like that on the level of 
equipartition in the PWN and that on how reliable is the 
age estimated from the pulsar spin-down properties. 



G54. 1+0.3 is the Galactic PWN that most c losely re- 
sembl es the Crab Nebula: this is the reason why iLu et al.l 
(|2002[ ) have dubbed i t "a close cousin of the Ciah Nebula" . 
Using Chandra data, iLu et all (|2002[ ) have shown the pres- 
ence of a well defined torus of ~ 10" in diameter, to- 
gether with elongations, toward E and W directions, which 
could be ascribed to X-ray jets. From those data, the size 
of the X-ray nebula appears 1', but the outer part of 
the nebula is very faint, and its edge is poorly defined. 
At radio wavelength s, instead, the nebular size is ^ 1.5' 
(jVelusamv &: Becker|[l98 8). corresponding to ^ 2.7^6.2 pc 
where de,2 is the distance of G 54. 1+0.3 in units of that 
estimated bv iLeahv et al.l (|2008[) . namely de,2 — d/6.2lj;g 
kpc. It is important to understand to which extent this dif- 
ference in size is real (i.e. due to synchrotron losses of the 
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Table 1. X-ray observations used in this work 



SateUite 


ID 


Texp (ks) Obs. date 


XMM-Newton 
SUZAKU 


0406730101 
502077010 


4r 26.09.2006 
8~ 30.10.2007 



PN exposure time after proton flare screening 



emitting electrons), or it is an artifact of the limited X-ray 
sensitivity. 

Radio maps show a rather amorphous structure, but 
the radio emis sion from G54. 1-1-0.3 is h ighly polarized, 
up to 20-30% ()Velusamv fc Beckeil [19881) . and this indi- 
cates (similarly to the case of the Crab Nebula) that the 
nebular field is highly ordered. The X-ray spectrum is a 
power law with a photon index —1.9, an absorption col- 
umn density Nh ^ 1.6 x 10^^ c m~^ and an X- r ay lum inos- 
ity Lx 1.3 X IQ^^dl ^ ergs"^ ICamilo et al.l (|2002|) have 
detected the pulsar PSR J1930-I-1852 at the center of the 
nebula, which has a period of 136 ms, a characteristic age 
of 2900 yr and a spin-down luminosity of 1.2 x 10"^^ erg s~^. 

Re cently, G5 4. 1-I-0.3 has attracted some interest be- 
cause iKoo et aT] (|2008l) have identified an IR shell sur- 
rounding the PWN at a distance of ~ 1.5' from the pul- 
sar. The shell c ontains a dozen of IR compact sources. 
iKoo et all (|2008f) suggests that the sources are young stel- 
lar objects, whose formation has been triggered by the wind 
of the progenitor of the SN. This i ntrigu ing possibility has 
been questioned bv iTemim et al.l (j2010t) . who pointed out 
that the IR shell may be ejecta dust, rathe r than a pre- 
existing ISM dense cloud. iLeahv et al.l (|2008[ ) reported the 
presence of a molecular cloud partially interacting with the 
PWN, on the basis of the CO emission around the nebula. 
Therefore, even if there is no hint of a radio shell around 
this PWN, there are a number of evidences for interac- 
tion between the PWN and the surroundings, so it is worth 
searching for an X-ray shell. 

In Sect. 2 we present a deep X-ray campaign aimed to 
this PWN, which led to the detection of such a shell. In 
Sect. 3 we estimate the contribution of the dust scattering 
halo, showing that it is negligible at the shell location, and 
in Sect. 4 we discuss our findings comparing them to a 
PWN-SNR evolution model. 

2. Deep XMM-Newton and SUZAKU observations 

We have observed t he PWN G54. 1+0.3 with XMM-Newton 
(jJansen et al.ll200l[ ) and (Mitsuda _ct^,2007i) in 

2006 and 2007 respectively. Table [T] summarizes the obser- 
vations used in this work. The data have been analyzed with 
the latest software available, namely SAS v8.0 for XMM- 
Newton and the pipehne v.2.1.6 with HEASOFT v6.6 for 
SUZAKU. The XMM-Newton data have been screened for 
pr oton flares using th e signi a clipping algorithm described 
in ISnowden k. Kuntd (|2007D . while for SUZAKU we have 
used the standard screening. The image of G54.H-0.3 ob- 
tained w ith the EPIC PN and MPS CCD camer as of XMM- 
Newton ^m^^^l^^^^^^^^^^M, with a res- 
olution of 8" HPD) a nd the XIS CCD camera of SUZAKU 
(jKovama et al.ll2007L 110" FWHM) are shown in Fig.[Il in 
the upper and lower panel respectively. 

There is a diffuse extended emission surrounding the 
bright plerion in both images. This diffuse emission seems 
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Fig. 1. Top panel: XMM-Newton EPIC image in the 1-7 
keV energy band. The image is a mosaic of background 
subtracted and vignetting corrected images obtained with 
the PN, MOSl and M0S2 cameras. The color scale has 
been chosen to maximize the visibility of the faint diffuse 
emission around the bright PWN. The inner (black) circle 
has a radius of 160" and includes the extension of the radio 
plerion, while the outer (white) circle has a radius of 400" 
and is drawn to guide the eye on the extension of the shell. 
Black thin contours marks an increase of factor 2 in surface 
brightness from 3x10^^ to 3.2 x lO^'^ cnt s~^ pix""'^, where 1 
pixel is 4" x 4". Red contours are drawn form t he map of the 
53 km s~^ CO emission published in Fig. 3b of lLeahv et al.l 
(|2008[ ). and represent a molecular cloud possibly interacting 
with the PWN. Bottom panel: SUZAKU image in the 1-7 
keV energy band. The image is a mosaic of background 
subtracted and vignetting corrected images obtained with 
the XISO, XIS2 and XIS3. The field of view and circles are 
the same as in the top panel. 



to have spatial structures as small as 1'. The contrast of 
surface brightness between the diffuse emission and the core 
of the PWN is about 1000. 

We investigated the nature of the X-ray faint diffuse 
emission by studying its spectrum both with XMM-Newton 
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Table 2. SUZAKU spectral fitting results 



Core of G54. 1+0.3 
power-law model. SUZAKU XIS+HXD 



Model 


N„ kT or 7^ 
10^2 cm-2 keV ~ 











Core of the PWN^ 








Power-law 


1.57 ±0.06 1.82 ±".04 


925/1145 


> 






Sheir 


CD 




Power-law 


1.57 "2.9 ± 0.3 


173/215 







Thermal 


1.57 2.0 ±0.4 


196/215 








" the power-law photon index 

Unabsorbed flux in the 2-10 keV band is 7.5 ± 0.5 x 
10"^^ erg cm"^ s"^ 

Extraction region 91 arcmin^. Unabsorbed flux in the 
2-10 keV band is 4.7 ±0.7 x 10"^^ erg cm"^ s"^ without 
the contribution scattered from the core. 



and SUZAKU. For the spectral analysis purposes, we de- 
fined a core region and a shell region. The core region is 
a circle centered on the pulsar with a radius of 160". This 
includes all the regions occupied by the radio nebula and 
its X-ray emission is dominated by the synchrotron emis- 
sion coming form the PWN. The shell region is an annulus 
with a inner radius equal to the radius of the core region 
and an outer radius of 400". We estimated that the con- 
tribution of the core emission in the shell region is 25% 
for SUZAKU and 1% for XMM-Newton, purely due to the 
instrumental PSF wings and excluding the additional con- 
tribution of dust scattered X-rays. In the case of XMM- 
Newton, the shell region was further reduced to a pie sector 
between polar angles 90° and 135° (from N, anti-clockwise, 
where most of the knots are located. The background was 
taken from the same chip for SUZAKU, while in the case 
of XMM-Newton we have used both blank fields and the 
same observation to collect background (in the latter case 
an annular region with Rmin = 400" and Rmax = 450" 
was used), verifying that the results did not change with 
the particular choice of the background. We fitted the core 
region, finding that a power-law model describes the data 
very well. Therefore, we fitted the shell region using a com- 
bination of the power-law model used in the core (with 
parameters fixed to their best-fit core values, including in- 
terstellar absorption) and an additional component, chosen 
among a thermal and a non-thermal model. The core model 
was used (and rescaled) in the shell region to take into 
account possible contamination from dust-scattering and 
from instrumental Point Spread Function. The SUZAKU 
results are summarized in Table [21 and seem to indicate 
that the shell emission which can be modeled either with 
a thermal component or with a non-thermal power law in 
the shell spectrum. XMM-Newton gives similar results. The 
SUZAKU spectrum of the core and the shell are reported 
in Fig. [2] along with their best-fit models. The core is also 
detected between 15 and 25 keV using the non-imaging 
SUZAKU Hard X-ray Detector (HXD) silicon PIN diodes 
(spectrum also reported in Fig.[51 upper panel), with a flux 
of ~ 4 X 10~^^ erg cm~^ s~^ in this band. The normaliza- 
tion constant between the PIN and the XIS CCD spectra 
is 1.3(1. — 1.6), a range which includes the expected value 
of 1.15 (iKokubun et al.l[2007lV 




2 5 
Energy (keV) 

Shell of G54.1+0.3 
Fixed (core) power-law + mekal. SUZAKU XIS 




Energy (keV) 

Fig. 2. Top: SUZAKU spectrum of the core region and the 
best-fit model. XISO+3, XIS2 and HXD spectra are shown. 
Bottom: SUZAKU XISO+3 and XIS2 spectra of the sheh 
region. The total best fit thermal + residual non-thermal 
component from the core is overplotted, as long as the in- 
dividual components (dashed). 



3. Removal of dust-scattering halo 

The presence of a foreground medium may affect in more 
ways the observed emission from an X-ray source. The best 
modeled effect is photoelectric absorption, but also dust 
scattering of the X-ray photons may be important. Its pro- 
duces an apparent halo around the intrinsic source, which is 
more prominent at lower energies and may hamper consid- 
erably both spectral mapping analysis of diffuse sources and 
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searches of faint surrounding features. Scattering hales may 
be relevant whenever the column density of the intervening 
material (Nh) is high. Nh may be derived by fitting the 
photoelectric absorption of the X-ray spectrum, while the 
scattering optical depth (rsca, estimate d at a reference pho- 
ton en ergy of 1 keV) by fitting the halo. lPredehl fc SchmittI 
([199^ show that the linear regression 

r,ca(l keV) = 0.05(7Vh/102i cm^^) ^ 0.083 (1) 

can be drawn between these two quantities. For the mea- 
sured value of 1.57 x lO^^cm"^ for G54.1-t-0.3, this relation 
predicts Tsca(lkeV) ~ 0.7. This value is however only ap- 
proximate, because the properties of the dust may be dif- 
ferent along different directions, while for the exact value a 
direct analysis of the halo is required. 

While halos of strong, point-like X-ray sources are se- 
lected to investigate properties of the dust grains, in this 
case we have to deal with a fainter, diffuse intrinsic source 
and therefore the modeling of the halo is necessarily less ac- 
curate. On the other hand, our goal here is simply to find a 
modeling that can subtract efficiently the halo component, 
without pretending to infer reliable physical properties of 
the dust distribution. 

In the case of a point-like source, we can describe the 
halo with the following function: 

IsciO, E) = F{E) (1 - exp(-r,ea(S)) H{e, E). (2) 

where F{E) is the source intrinsic flux and H{0, E) is 
a function that we derive and describe in detail in the 
Appendix]^ where we also summarize the basic modeling 
of a halo and the related assumptions. 

In our case, however, the intrinsic source is spatially re- 
solved, and in principle we should convolve the halo for a 
point-like source with the surface brightness distribution of 
the actual intrinsic source, in a way to obtain the observed 
map. This is in general a very complex and numerically 
heavy task. Moreover, if just one energy band is consid- 
ered, the number of possible solutions would even be in- 
finite. To solve the problem of separating intrinsic source 
and halo, one must then carry on a combined fit on radial 
profiles at different energies, by taking advantage of the 
known energy dependence of the halo properties. In order 
to model the energy dependence of the scattering halo, we 
have produced for each instrument 4 images, in the follow- 
ing spectral bands: 1.0-1.75 keV, 1.75-2.5 keV, 2.5-3.9 keV, 
and 3.9-7.25 keV. The spectral boundaries have been cho- 
sen to get similar numbers of photons in the various bands, 
with the further constraint of excluding photons softer than 
1 keV, for which the simple Rayleigh-Gans scalings with the 
photon energy (namely optical depth oc E^"^ and halo size 
oc E~^) are no longer valid. The reference energies of the 
four bands are 1.4, 2.05, 3.07, and 4.94 keV respectively, 
and are obtained averaging the energies of all photons col- 
lected in each band. 

Since the intrinsic source is centrally peaked and more 
concentrated than the halo, and the wings of XMM-Newton 
PSF are narrower than the observed radial profile, we have 
applied a simplified approach, by approximating the halo 
with that for a point-like source. For our final fits we have 
used only 2 free parameters to model the halo, namely 
rsca(l keV) and 6'scai(l keV). Before then, we had attempted 
also fits using a larger number of parameters, but with the 
moderate statistics of our data we have found: i. a par- 
tial degeneracy between power-law index of the grain size 



Table 3. X-ray dust scattering halo best-fit parameters for 
054.14-0.3 



Name 


0" - 120" 


0" - 200" 




6.8" 


6.6" 


A 


0.87 


0.97 




30.7" 


29.2" 


0-2 


2.27 


2.22 


T,ca(lkeV) 


1.15 


1.09 


6i,cai(lkeV) 


15.1' 


10.2' 



distribution (g) and the spatial scale (^scai), and therefore 
we have chosen a rathe r usual value for q (3.5, see e.g. 
iPredehl fc Schmittlll995D : ii. our fits were typically consis- 
tent with a wide spread of z (i.e. the position along the line 
of sight, normalized to the source distance), and therefore 
we have decided to assume a homogeneous distribution of 
dust along the line of sight (zmin = and z^ax = !)• AH 
quantities cited here are described in the Appendix |^ 

For the intrinsic source we have chosen the following 
modeling: 

W^, i^) = ^ (exp (-|) + A (1 + I) , (3) 

where C = 2al + Aa"^/ {'^{0.2 — 1)) is the normalization fac- 
tor. The fitted shape is indeed a convolution of the actual 
source with the instrumental PSF, but at this level we do 
not need to separate theirQ- In addition, we have assumed 
that this shape is independent oi E: indeed, it is known that 
the X-ray size of PWNe is slightly decreasing for increasing 
E, but this is only a minor effect, which we cannot ade- 
quately describe with the available data, and on the other 
hand does not affect considerably the results of our fits. 

From each image we have extracted a logarithmically 
spaced radial profile. In order to minimize the statistical 
noise, the flux values are averaged over several points. In the 
case of MOSl, due to the absence of the damaged CCD#6, 
data are missing for a region that is relevant to our pur- 
poses, so we use M0S2 only. By analyzing the emission at 
distances larger than 300" from the source center, we have 
estimated the M0S2 background levels for the 4 bands, as 
about 10""' times the surface brightness in the brightest 
areas. Therefore, we did not apply any correction of the 
background. With this, we are confident that our profiles 
are usable over a dynamical range close to 10^. 

For the fits, we have used at first only the inner 120" 
of the radial profiles. This to allow an unbiased analysis of 
structures that may appear at larger radii. However, we 
also present the results obtained in the inner 200" of the 
radial profiles. The results of the fits are presented in Table 
[21 The results are in general rather similar between different 
choice of maximum fitting radius. 

Fig. [3] shows, the combined M0S2 fits on the 4 bands 
(for a maximum fitting radius of 120" marked by the verti- 
cal dashed line) and the extrapolation of the best-fit pro- 

^ Eq. 3 corresponds to the PSF analytical descrip- 
tion discussed in the EPIC Calibration status document 
available in the ESA XMM-Newton Cahbration Portal 
(http://xmm2.esac.esa.int). However, our best-fit parameters 
are in general larger since the PWN is extended. To the pur- 
pose of our halo modeling, we just need an empirical relation to 
take into account the PSF+source effect. 
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files to larger radial distances, compared to the observed 
profiles. The observed radial profile shows an excess of 
emission at distances from the center > fOO" which cannot 
be explained by the dust-scattering halo. If we normalize to 
the fluxes due to the source plus the halo in the 0" — 160" 
region, the model predicts that the flux due to the halo 
in the outer 160" - 400" annulus are 0.20, 0.10, 0.04 and 
0.01 in the 1.0-1.75, 1.75-2.5, 2.5-3.9 and 3.9-7.25 keV re- 
spectively (the flux of the source is 0.01 in the same re- 
gion), while the observed values are 0.28±0.01, 0.14±0.02, 
0.11 ± 0.02, 0.16 ± 0.03, significantly higher than the pre- 
dicted halo model. We therefore conclude that the excess is 
due to intrinsic emission from the G54. 1-1-0.3 shell. 

4. Discussion 

We have reported for the first time the detection of a large 
faint diffuse X-ray emission around the PWN G54.1-I-0.3. 
We have seen that this excess cannot be due to the X- 
ray dust scattering halo, because an accurate modeling of 
the halo presented in previous section and in the Appendix 
A shows that the halo model underpredicts the observed 
emission between 160" and 400" from the center, where the 
shell is observed. This shell emission has an irregular mor- 
phology, but can be enclosed at most inside a circle of 5.7 
arcmin radius centered on the pulsar position, correspond- 
ing to ^ 10.3 C?6.2 pc. 

We have seen that the halo spectrum can be inter- 
preted as thermal or non-thermal emission. If we interpret 
it as the long sought thermal emission of the G54. 1-1-0.3 
shell, we can derive some interesting quantities related to 
the re mnant evoluti on by assuming an expansion governed 
by the ISedovl (|1959D solution. The best-fit emission mea- 
sure is EM = 9 X 10~^^ cm~^, and the best-fit tempera- 
ture is 2 keV. According to the model of iGhavamian et al.l 
(|2007( ). electrons cannot be heated directly at this tem- 
perature by the shock, so the Coulomb heating by the 
shocked ions must be at work. Using the relations (6)-(10) 
of Bocchino & Bandiera ( 2003), and the best-fit emission 
measure and temperature, we have computed the plot of 
Fig. m which shows all the possible solutions versus the 
remnant distance. We have specifically taken into account 
the case in which electrons are not thermalized with ions, 
by reporting 3 different cases of electron-ion temperatures 
in Fig. 01 namely Te/Tp = 1 (electron-ion full equilibra- 
tion), Te/Tp = 1/2 (moderate electron- ion disequilibrium) 
and Tf lTr, = 1/10 (strong dise quilibrium, as in other young 
SNRs. lGhavamian et al.|[2007l and references therein). If we 
use a distance of 6.2 kpc, we find a range of remnant ages 
between 2500 and 3300 yr for the equipartition case (solid 
line in Fig. U]), which is in agreement with the more uncer- 
tain estimate of the remnant age (1500-6000 yr) given by 
ICamilo et all (|2002[ ). The inferred ISM pre-shock density 
is ~ 0.2 cm~^ and the swept up mass is between 23 and 
32 Mq, while the X-ray emitting mass is 15-20 Mq. The 
explosion energy range is i? = 3 — 7 x 10^" erg. However, if 
we drop the assumption of equipartition between electron 
and ions, we find that E = 1.0(0.5 — 1.6) x 10^-^ erg for 
a distance of 6.2 kpc when T^/Tp ~ 1/2. In this case, the 
derived age is between 1800 and 2400 yr (dotted line in Fig. 
[3]). Lower values of Tg/Tp (i.e. ^ 1/10) are disfavored by 
relative high explosion energies (dashed line in Fig.|4]). We 
have seen that the XMM spectral fittings suggest a similar 
temperature to the SUZAKU values, but a normalization 



10 times lower. In this case, the estimate of the ISM den- 
sity, explosion energy and swept-up mass must be decreases 
by a fact or of 3. A cross-chec k with the non-radiative SNR 
model of iTruelove fc McKed ([1999) gives a transition from 
ejecta-dominated to Sedov phase at 2500 yr, so the rem- 
nant is entering in the adiabatic phase. Given the faintness 
of the diffuse emission, only a deeper X-ray observation 
would allow us to derive more reliable values of the shell 
parameters. 

The derived values of density and swept-up mass suggest 
an expansion in rarefied medium for most of the remnant 
lif etime. This s eems to be in agreement with the findings 
of iLeahv et al.l (|2008D about the environment of G54. 1-1-0.3. 
The remnant projected location inside a large IR shell opens 
up the possibility that G54. 1-1-0.3 have been originated by 
a SN belonging to the same star cluster whose winds have 
created the large IR shell. Although the IR shell distance 
seems to be a little bit large r than the PWN di stance (7.2 
kpc vs 6.2 kpc), according to lLeahv et al.l (|2008l ) the uncer- 
tainties on the distance do not rule out the association, and 
the values of the density we derived with the X-ray spectral 
analysis goes in the same direction. The X-ra y shell seems 
to be larger than the CO cloud reported by iLeahv et al.l 
1 2 0^81). as shown b y the CO contours overplotted in Fig. 
□ Koo et all (|2008 h showed that G54.1-^0.3 has interacted 
with a star-forming loop located very close to the nebula 
center (at ~ 1' from the pulsar). The loop contains at least 

11 young stellar objects (YSOs) which are very b right in the 
AKAR I 15 /Ltm image of the core of the PWN. iKoo et al.l 
(|2008D argue that there is no direct evidence for the inter- 
action of the SNR shock with this dense material, so they 
conclude that the IR loops is a partial shell in a low-density 
medium and that the SNR shock has propagated well be- 
yond it. This is in agreement with the position of the X-ray 
shell w e have discovered, since we can now compute (us- 
ing the ITruelove fc McKe3ll999l model) that the shock was 
at the IR lo o p pos ition just at 1/10 of the present age. 
iTemim et al.l (|2010l ) proposed an explanation in terms of 
ejecta dust for the IR loop, which is not in contrast with 
the presence of the X-ray shell. 

If the faint diffuse emission we have discovered around 
the PWN G54.1-h0.3 is really the associated SNR shell, 
then we can compare its pro perties with th e PWN-SNR 
evolutionary model of Gelfa nd et al.l (j2009t ). This model 
couples the dynamical and radiative evolution of the PWN 
with the dynamical properties of the surrounding non- 
radiative SNR, and predicts several distinctive evolution- 
ary stages, namely the initial expansion, the reverse shock 
collision, the re-expansion and the second compression. We 
run the model using as input parameters i?sn, 51 = 1 and 
nism = 0.2 cm~^, derived from the best-fit thermal model 
of the X-ray sheU (Fig. HI T^/Tp = 1/2 case of a moder- 
ate deviation from electron-ion equipartition) . We assumed 
an ejecta mass of 8 Mq, and a spin-down timescale of 500 
years (close to the value of the Crab; this corresponds to 
an initial period Pq = 56 ms of the pulsar). The pulsar 
wind p roperties are the same as in Table 2 of Gelfand et al] 
(I2OO9I ). The resulting dynamical evolution of the shock of 
the shell and the PWN nebula is shown in Fig. |5j where 
we plot the ratio between the shell radius and the PWN 
radius {Rsnr / Rpwn) versus time, and a vertical line marks 
the time when t he pulsar has the sa me characteristic age 
as measured by ICamilo et al.l (|2002l) , and its period and 
period-derivative match what was observed (tc = 2900 yr 
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and a real age of 2400 yr). The predicted Rsnr / Rpwn is re- 
markably similar to the observed value, and the SNR and 
PWN sizes (8.7 and 2.8 pc) are in relatively good agreement 
with observation^ According to this model, the PWN has 
not yet been crushed by the reverse shock (Fig.[5]shows that 
it will happen at an age of ^ 7000 yr, when Rsnr / Rpwn 
reaches a minimum), in agreement w ith the lack of sign 
of crushing, as noted bv lTemim et al.l ([2010). We conclude 
that the thermal parameters we have measured in the faint 
diffuse emission around the PWN and its size are in good 
agreement with what we expect from a putative shell of the 
SNR G54. 1-1-0.3, on the basis of a complete modeling of the 
SNR-PWN system. We interpret this as a further indica- 
tion that what we observe is indeed the long sought shell 
of the remnant. 



5. Summary and conclusions 

We have analyzed an XMM-Newton and a SUZAKU obser- 
vation of the PWN G54. 1-1-0.3, in the framework of a pro- 
gram aimed to survey the region around this isolated neb- 
ula in search for the X-ray shell of the associated supernova 
remnant. We detected very faint X-ray emission around the 
PWN, extending form the outskirts of the PWN (at 1.5' 
from the central pulsar) until a radius ^3.8 times the PWN 
radius (i.e. ~ 5.7', around 10.3 pc at the distance of the neb- 
ula) . This extended diffuse emission is more evident toward 
south and it has an irregular morphology on a angular scale 
of ~ 1'. We modeled the X-ray dust scattering halo around 
G54.1-(-0.3, and we have found that the detected faint dif- 
fuse emission cannot be due to this effect, but it must be 
intrinsic to the source. We modeled the X-ray spectrum of 
the diffuse emission with a thermal model, finding a best-fit 
temperature of ~ 2 keV, which may imply electron heating 
by the shocked ions. This value, together with the appar- 
ent size and the emission measure of the X-ray emitting 
plasma, is consistent with a SNR shell expanding into a 
~ 0.2 cm~'^ ISM, whose explosion energy is ~ 10^^ erg, 
and whose most probable age is 1800-2400 yr, a bit less 
than the characteristic age of the pulsar PSR J1930-I-1852, 
located at the center of the PWN. However, due to limited 
counting statistics, the X-ray spectrum of the diffuse emis- 
sion can be alternatively well fitted with a non-thermal 
power-law model, whose photon index (7 = 2.9) is roughly 
consistent with an interpretation in terms of synchrotron 
emission from accelerated particles. The morphology of the 
large diffuse emission neither seems to be dire ctly linked 
to the IR shell observed around the PWN b v iKoo et al 



( 2008 ). nor to the molecular cloud detected bv lLeahv et al 



2008[) and reported as contours in Fig. [U but the fact that 
the X-ray shell is incomplete is probably related to the in- 
teraction between the PWN and these inhomogeneities of 
the ISM. 

We have compared the PWN and SNR siz es with 
the pr ediction of the evolutionary model of Gelfa nd et al.l 
(|2009t) for composite SNRs, and we find an excellent agree- 
ment. We conclude that the faint diffuse emission around 
the PWN G54. 1-1-0.3 may indeed be the shell of the associ- 
ated remnant. However, deeper X-ray and radio observa- 



We have verified that, by running the PWN-SNR model us- 
ing the results of the equipartition case (Te/Tp = 1 in Fig.[4|) we 
would not reproduce the observed size of the remnant, unless a 
far too low ejecta mass is used. 



tions are required to definitely distinguish between thermal 
and non-thermal interpretation. Given the recent detec- 
tions of X-ray shells around other PWNe, our results sug- 
gest that the lack of shell around remaining isolated PWNe 
may be simply the result of high absorption and/or lack of 
long observations, and that the X-ray band may be very 
effective in discovering them. 
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Appendix A: An approximated treatment of 
dust-scattering X-ray halos 

We assume that scatterings are in Rayleigh-Gans 
regime (this i s typi cally valid above 1 keV; see, e.g., 
ISmith fc Dwe^ (|19^a^ for more details on the different scat- 
tering regimes). In addition, we assume that scattering 
angles are small, and we do not consider multiple scat- 
terings fe.g.. iPredehl fc Klosd [19961) . In this way, one can 
take advantage of two simple scaling laws: the scattering 
optical depth scales with the photon energy as E~'^, while 
the angular scale of the halo scales as E~^. 

The radial profile of the halo of a point-like source, at a 
given photon energy, is derived by calculating the following 
integral 

IUO,E) = F{E)NhJ ^f{x) J dan{a)^{E,a,^), (A.l) 

where 9 is the (angular) radial distance, F{E) is the source 
intrinsic fiux, x is the normalized distance to the sourc e 
(x = 1 — z, following the notation of lSmith fc Dwel3ll998[ ). 
f{x) is the normalized distribution of dust density along 
the line of sight, and n{a) is the normalized distribution 
of grain sizes (a; here we assume a position-independent 
shape of this distribution). Finally, the scattering angle 4> 
is equal to 6/x. 

The differential cross section for single scattering can be 
expressed as 



^{E,a,^) = Ca'G{yr, 



(A.2) 



where C contains the information on the grain composition, 
y is defined as aE(f)/hc, and 



G(y) = ^(My) - ycos(y)) = 1 - ^ + 0(2/")- (A.3) 
y-^ 10 



Therefore, assuming a constant f{x) between Xmin and 
a^max and zero elsewhere, and a power-law grain size dis- 
tribution (n(a) oc a^"^ up to a maximum size amax, we have 

In the limit Omax ^ hcx/E6, the integral in a is equal to: 



while, in the opposite limit, it is equal to: 



^ hex 
EO 



7-q 



^scal 



7-9 



(A.5) 



(A.6) 
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where we have defined ^scai = C^c/amax-E and 



rf2/(7-g)y'-'G(y) 



(A.7) 



The quantity ^ is a function of q only: it evaluates 2.418 
for q = 4, 2.727 for q = 3.5, while it slowly diverges for q 
approaching 3. Ju st as an orientative va lue, for q — 3.5 and 
Omax — 0.25 //m (jSmith fc Dwekl [l998D we have Sgcai — 
lAE^Jy arcmin. 

We introduce a "step-like" approximation for the func- 
tion G{y)^: namely, equal to unity for y < ^ and vanishing 
elsewhere. This approximation is equivalent to approximate 
the integral in a by matching its two limits. In this way, it 
is possible to integrate analytically the integral in Eq. A.7. 
The result is proportional to function H{6.E) whose shape 
defined as 



(6 

(7- 

{x 



-9)(^min 

-(6 - g)a;„,^^ 



max/ 
''mill ''n 



if 6I„ < Xn 



6-q 



„S-q\ 



-(7"-3) 



if Xniin 

if Xry^a-V 



< ( 

< 



(A.8) 



where 9^ — d/9sca\- The total intensity of the profile as 
defined by Eq. A.7 is 



^ ~ {5-q) ^ " 



2 

seal" 



(A.9) 



It is then convenient to redefine H{6) as divided by W, so 
to normalize its integral. In this way, we can finally simply 
write: 
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I E) = F[E) (1 - exp(Tsca(£;)) H{e, E). 



(A.IO) 



With respect to the profile derived here, an exact solution 
would give a slightly smoother profile. However, this ef- 
fect is minor compared to the uncertainties related to how 
sharp is the upper cutoff in the distribution with size and 
how sharp are the boundaries of the spatial distribution of 
grains. 

It is worth noticing a few properties of H{9^ E). As ex- 
pected, 0sciii oc E~'^; also ^scai oc a^l^. In addition, if Xmin 
and Xmux are both changed by a factor 77, the shape of H 
does not change, provided that ^scai changes as 77. This, 
for instance, implies that the effect of dust extending from 
the source to a minimum normalized distance (zmin) to the 
observer is equivalent to the case of a uniform spatial dis- 
tribution of the dust, but with a size distribution extending 
to amax/^min instead of to flinax: this effect would lead to 
overestimate the value of flmax- 

Since here we do not want to study the actual properties 
of the foreground dust but simply to model the shape of the 
scattering halo, without loss of generality, in the following 
we assume — (i.e., Xmax = !)■ 
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Fig. 3. X-ray normalized profiles in the bands 1.0-1.75, 
1.75-2.5, 2.5-3.9, and 3.9-7.25 keV of G54.1+0.3 as seen 
by M0S2 camera of XMM-Newton EPIC. In the x-axis we 
report the logarithmic distance from the center in arcsec. 
We have overplotted the best fit halo model, with indi- 
vidual components (intrinsic source, long-dashed line, and 
halo, short-dashed line). The radial range of the fit extends 
only to 120" (the position of the vertical dashed line). See 
the text for further details on the model. The residuals be- 
tween 160" and 400" (2.2 to 2.6 in logarithmic scale) with 
respect to the extrapolation of the halo model are due to 
the X-ray shell. 
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Fig. 4. Sedov solutions for the SNR age and distance com- 
puted using the X-ray derived parameters of G54.1+0.3 
{kT = 2 keV, SNR radius 6', emission measure 9 x 10^-'^ 
cm~^). We report 3 different solutions corresponding to 
full electron-ion equipartition (Te/Tp = 1) and the cases 
oiTe/Tp = 1/2 and 1/10. The asterisks on the curve mark 
the corresponding explosion energy of the solutions. The in- 
tersection of the curves with the area marked in green (cor- 
responding to the distance estimate of lLeahv et al. (.2008) ) 
gives the age of the system (for instance, the age is 1800- 
2400 yr for the case Te/Tp = 1/2, and the energy is ~ 10^^ 
erg). 



Evolution of the PWN/SNR G54.1+0.3 
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Fig. 5. Ratio between the shell radius and PW N radius 
according to the model of iGelfand et al.l (|2009[ ). with in- 
put parameters tailored to the G54. 1-1-0.3 system (see text 
for details). The observed value is also reported as a cyan 
stripe, whose spread is due to the uncertainties in defin- 
ing the real PWN radius (we have used a minimum of 80" 
and a maximum of 114", on the base of the morphology 
shown in Fig. [IJ. The vertical dot-dashed line marks the 
instant of the run when the simulated pulsar has a char- 
acteristic age and spin-down luminosity identical to PSR 
J1930-hl852 as measured bv ICamilo et all (|2002[ ). In this 
moment, the model reproduces both the PWN and SNR 
size and their ratio, strengthening the conclusion that the 
extended diffuse emission is indeed the SNR shell of the 
PWN G54.1-I-0.3. 



